Spectrum of Kinetic Alfven Turbulence 



Stanislav Boldyrev 1 and Jean Carlos Perez 1,2 
1 Department of Physics, University of Wisconsin, Madison, WI 53706, USA 
2 Department of Physics and Space Science Center, U. New Hampshire, Durham, NH 
(Dated: 9 February 2012; revised 10 April 2012) 



USA 



(N 

o 

(N 

< 



9 L,: 
6 

CO 



> 
o 

00 

in 
o 

(N 



A model for strong kinetic Alfven plasma turbulence at scales smaller than the ion gyroscale 
is proposed. It is argued that magnetic and density fluctuations are concentrated mostly at two- 
dimensional structures, which leads to their Fourier energy spectra E(k±) <x fc^ 8 ^ 3 , where k± is 
the wave-vector component normal to the strong background magnetic field. The results are shown 
to be in good agreement with numerical simulations, and they can explain recent observations of 
magnetic and density fluctuations in the solar wind at sub-proton scales. 



PACS numbers: 52.35.Ra, 95.30.Qd, 96.50.Tf, 96.50.Ci 



Introduction — Possibly the simplest description of 
magnetic plasma turbulence at scales much larger than 
typical micro-scales (particle gyro-radii, skin depth, etc) 
is provided by one- fluid magnetohydrodynamics (MHD), 
where electrons and ions are assumed to move together 
as a single fluid [e.g., [H 0. A characteristic feature of 
MHD turbulence is its anisotropic spectral energy trans- 
fer with respect to the background magnetic field. As a 
result, small-scale plasma fluctuations populate predomi- 
nantly field-perpendicular wavevectors and turbulence is 
dominated by the shear- Alfven modes Q. In the linear 
case these modes have the dispersion relation uj — k^VA, 
where ku is the field-parallel wavenumber with respect to 
the background magnetic field Bo, va = Bo/^/4¥p is the 
Alfven velocity, and p is the fluid density. 

At scales smaller than the so-called kinetic- Alfven dis- 
persion scale (ion acoustic radius if the electron tempera- 
ture exceeds the ion temperature, or ion gyroradius oth- 
erwise) the assumptions of one-fluid MHD breaks down 
and the nature of turbulence changes. At such sub- 
proton scales, the shear-Alfven cascade transforms into 
the cascade of strongly anisotropic kinetic- Alfven modes 
with a different linearized dispersion relation ui oc k\\k±. 
Kinetic Alfven turbulence attracts considerable interest 
due to its importance for solar wind heating, magnetic 
reconnection in a variety of astrophysical systems, and 
laboratory experiments with strongly magnetized plas- 
mas [e.g.. I34l0|. Such turbulence has been understood 
to a much lesser extent compared to MHD turbulence. 

In this contribution we address the spectrum and struc- 
ture of kinetic Alfven turbulence. Our consideration is, 
in part, motivated by measurements of small-scale fluc- 
tuations in the solar wind. Although significant scatter 
exists among the reported data 



111, recent observations 



suggest that magnetic fluctuations at sub-proton scales 
have the Fourier energy spectrum close to or possibly 
steeper than fc -2 8 [e.g.. Hol. [r3 - [l5j . The nature of such 
fluctuations is unclear as they are not described by ex- 
isting models of either kinetic-Alfven or electron mag- 
netohydrodynamic turbulence, which predict the scaling 
fc~ 7 / 3 . Recently proposed explanations include signifi- 
cant steepening of the spectrum due to Landau damp- 
ing, presence of weak turbulence, wave-particle scatter- 



ing, etc. [e.g., 1M LL71. 

In this paper we analyze small-scale kinetic-Alfven 
turbulence using a two-fluid plasma description, which 
by its nature, does not take into account Landau 
damping and other wave-particle interactions. We found 
that the steeper than fc~ 7 / 3 energy spectrum persists in 
this case, in the form closely resembling the solar wind 
observations and the results of existing kinetic simula- 
tions. Our results indicate that the power-law energy 
spectrum of strong kinetic-Alfven turbulence in not an 
artifact of significant dissipation or non-universality, 
but rather an inherent property of nonlinear plasma 
dynamics. To describe this spectrum we propose a 
new model that assumes that the magnetic and density 
fluctuations tend to spontaneously organize into two- 
dimensional structures thus leading to strong spatial 
and temporal intermittency of turbulence. This hap- 
pens as a consequence of kinetic-Alfven dynamics that 
combines nonlinear striation in the field-perpendicular 
direction and linear spreading in the field-parallel direc- 
tion. Our model predicts that the energy spectrum of 
strong kinetic-Alfven turbulence has the Fourier scaling 
oc k^ , which is in good agreement with our 
numerical findings and may explain the solar wind data. 

Kinetic Alfven model. — Let us assume that a uniform 
background magnetic field (the guide field) is strong com- 
pared to magnetic fluctuations, and we are interested in 
frequencies smaller than the ion gyrofrequency. Since 
the electrons are strongly magnetized and their thermal 
speed exceeds the Alfven speed, an isothermal fluid de- 
scription is possible for the electrons [Uj]. The electrons 
are advected across the guide field by the "E cross B" 
drift, vi = cE x B0/-B 2 . while their field-parallel motion 
is related to the current J|| = — en e i> e ||, and the ion par- 
allel motion can be neglected. For simplicity, we start 
with the case when the thermal plasma energy is small 
compared to the magnetic energy (small plasma "beta"). 
In this case the field-parallel fluctuations of the mag- 
netic field can be neglected, while its field-perpendicular 
component is expressed through the flux function = 
z x Vip, so that J|| = (c/4n)V± xb ± = (c/4tt)V 2 l -0. The 
flux function is the field-parallel component of the vector 
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potential, ip = —A z . 

The field-parallel force balance in the electron momen- 
tum equation gives — r e V||7i e — noeE|| = 0, where the 
electric field is E = — — (l/c)djA. Supplementing 
this equation with the electron continuity equation, one 
obtains the system for the fluctuating parts of magnetic 
and density fields: 

l^-V||^+^V||n e = 0, (1) 
| ne _^ V 0x|.Vn e -|V||J||=O, (2) 

These equations have been derived and studied in many 
works, e.g., [T§-23]. The field-parallel gradient in these 
equations is the gradient along the local magnetic field, 
that is, 

V|, = V 2 + ^x W-V. (3) 

-DO 

We will assume that the fluctuations are anisotropic with 
respect to the local magnetic field, so that the so-called 
critical balance between the linear and nonlinear terms 
is satisfied, V z ~ (1/B )z X Vip • V, [e.g., 0, i, E3 Hi- 
This is the case of strong turbulence that we consider 
in this paper. Equations ([TJ [5]) are therefore essentially 
nonlinear and three-dimensional. 

We still need to specify the electric potential <p in the 
system ([TJ [2J . We are interested in the dispersive kinetic- 
Alfven waves, that is, we consider the scales smaller than 
the ion- acoustic scale k±p s ^> 1 or the ion gyroscale 
if T e ~ Ti] here p s = C s /Q h C s = (Te/rm) 1 / 2 is the 
ion acoustic speed, and fli is the ion gyrofrequency. Be- 
low the ion gyroscale the ions are not magnetized, and 
since we are interested in frequencies smaller than kvxi 
(vn = {Ti/ mi) 1 / 2 is the ion thermal speed), we have 
for the ion density fluctuations ni = —ecpno/Ti. In this 
case, the quasi-neutrality condition rii = n e ensures that 
the second (advection) term in Eq. ([2]) vanishes, while in 
Eq. (p} the electric potential modifies the density term: 
V||0 = -(T i /n o e)V||n e HI. 

Let us introduce the normalized electron density 
h = (1 + Ti/T e ) 1 / 2 {C s /vA}n e /nQ^ magnetic flux func- 
tion ip = (C s /c)eA z /T e , and the electric potential <p = 
(C s /vA)e<fi/T e . We normalize the spatial scales to the 
ion-acoustic scale p s , and the time scale to (p s /va)(1 + 
Ti/Te)- 1 / 2 . In what follows we will use only the nor- 
malized variables and omit the over-tilde sign. Then we 
obtain that the magnetic and density fields in the kinetic 
Alfven regime are described by the system: 

ft^ + V||n = 0, (4) 
d t n - V||ViV = 0, (5) 

where V|| = W z +zxWtp-W ±. The presented ideal system 
conserves the total energy E and the cross-correlation H , 

E = J (|V?A| 2 + n 2 )d 3 x, (6) 
H = J i>nd z x. (7) 



The system (|4I5[) possesses linear waves, nt oc ipk oc 
cxp(— iut + ikx). The linearization is done by neglecting 
the second term in the right-hand side of Eq. (J3J , which 
gives the dispersion relation for the kinetic Alfven waves: 

uj = /c||fc_L. (8) 

The linear modes are characterized by the equipartition 
of density and magnetic fluctuations, — ±kj_ipk- 

To conclude this section we make two important com- 
ments. First, a similar consideration can be conducted 
without the assumption of small plasma beta. In this 
case the field-parallel fluctuations of the magnetic field 
should be taken into account in the derivation of ([TJ 
12) [ e -g-> US IHl- The resulting system however has the 
structure identical to our system ((H [5]) and it can be 
reduced to system (|4j [5]) by appropriate normalization 
of the variables. The value of beta is therefore not es- 
sential for our discussion of scaling properties of kinetic 
Alfven turbulence. Second, we are interested in the 
processes slower than the ion gyrofrequency, implying 
fey 1. However, equations ((H [SJ admit a rescaling 
d/dt — > ed/dt, V 2 — > eV 2 , n — > en, and tp — > op 
with arbitrary e, which preserves the critical balance. 
This reflects the fact that being derived in the limit of 
infinitely large electron gyrofrequency Eqs. (jH [5]) lack 
any frequency scale. We may therefore always rescale the 
fields in these equations to satisfy fen ~ n ~ tp ~ 1. Such 
rescaling will be used in our numerical simulations below. 

Kinetic Alfven turbulence. — The scaling of strong ki- 
netic Alfven turbulence was addressed in a number of 
works, [e.g., HI, HI], see also 0, H H3, [H| . It was argued 
that in strong turbulence, the critical balance condition 
which ensures that both linear and nonlinear terms in Q 
are of the same order, should be satisfied at all scales. 
Denote n\ and ip\ the typical (rms) fluctuations at the 
field-perpendicular scale A, and I the corresponding field- 
parallel scale of those fluctuations. Balancing linear and 
nonlinear terms in ([3]) then gives I ~ X 2 /ip\, in which case 
the time of nonlinear interaction is comparable to the 
linear time ©, t ~ 1/w ~ ZA ~ X 3 /tp\. Besides, we esti- 
mate from that n\ ~ The energy associated 
with the scale A can therefore be estimated as E\ ~ n\, 
and the condition of constant energy flux in the turbulent 
cascade leads to n\/r = const, which translates into the 
scaling for the turbulent fields n\ ~ ip\/X ~ A 2 / 3 . The 
Fourier energy spectrum of strong kinetic Alfven turbu- 
lence is then: 

E KA {k ± )dk ± oc fc~ 7/3 <ife ± . (9) 

As we discussed in the introduction, there is a puzzling 
disagreement of this scaling with the solar wind observa- 
tions, where a spectrum closer to —2.8 is observed. 

To address this issue we have conducted numerical 
simulations of system (|4I5[) . Our results produce a 
turbulent spectrum that is different from ([5]), and 
quite close to the observational data. Since our system 
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FIG. 1: The energy spectrum of strong kinetic Alfven turbu- 
lence below ion acoustic scale (ion gyroscale). 

does not include Landau damping and it is driven in 
the regime of strong turbulence, we propose that the 
observed scaling is not an artifact of non-universal or 
dissipative effects, rather, it is an inherent property of 
nonlinear turbulent dynamics. We then propose a model 
of kinetic Alfven turbulence, which predicts that the 
energy spectrum should scale as E(k±) oc k ± 8 ^ 3 , in good 
agreement with our numerical results. 

Numerical simulations of kinetic Alfven turbulence — 
We supplement the system (|4I5[) by a driving force and 
by small dissipation terms as follows: 

d t ^ + V\\n = nV\^ + f, (10) 
d t n - V\\V 2 x ip = vV\n. (11) 

The force mimics energy supply from large-scale motion, 
while the dissipation terms (normalized plasma resistiv- 
ity r\ and electron diffusivity v) remove the energy at 
small scales; the dissipation terms are also needed to en- 
sure numerical stability of the code. We solve these equa- 
tions on a triply periodic cubic domain (L 3 , L = 1) using 
standard pseudo-spectral methods. The random force / 
is applied in Fourier space at wavenumbers 2tt / L < fcj_ < 
2(2ir/L), fc|| = 2tt/L. The Fourier coefficients outside the 
above range are zero and inside that range are Gaus- 
sian random numbers with amplitudes chosen so that 
IVV'Irms ~ 1- The individual random values are refreshed 
independently on average every r = 0.1L/(27r| V"0| rTO5 ). 
We choose v = rj = 0.01. The strength of the nonlinear 
term relative to the dissipation term is then measured by 
the parameter R = ip rms /(Li>) ~ n rms /v, which plays a 
role of the Reynolds number in this system. 

We use numerical resolution of 512 3 collocation points. 
The initial conditions are imported from a steady state 
snapshot obtained on 256 3 points. The system is then 
evolved until a new steady state is reached. The 512 3 
simulations are run for about 35 large-scale dynamical 
times. The presented results correspond to statistical 
averages over approximately 60 last snapshots cor- 
responding to about 15 dynamical times. Note that 
compared to the MHD equations where ui ~ fcii , the 
kinetic Alfven equations require significantly shorter 



time steps to accommodate high frequencies (|SJ), leading 
to tremendous increase in computational effort. In this 
respect the fluid model p0lll[) allows one to access 
the inertial intervals and averaging times currently un- 
achievable in kinetic or gyrokinetic simulations. Fig.([T]) 
shows the energy spectrum of kinetic Alfven turbulence. 
The spectrum is steeper than —7/3 and close to —8/3. 

A model for kinetic Alfven turbulence — To understand 
the observed energy spectrum, let us discuss some char- 
acteristic properties of the dynamics described by the 
system ([HE]). First, consider the effect of the nonlinear 
terms (for that we can assume k\\ — 0). The nonlinear 
term in Eq. (|10|) can be rewritten as Vn x z ■ Wtp, imply- 
ing that i\) is advected in the field-perpendicular direction 
with velocity Vn x z. The field ip thus gets striated, de- 
veloping gradients aligned with the gradients of n. This 
suggests that the magnetic field (Vi/0 tends to concen- 
trate in 2D structures. Let us now see what happens to 
those structures if the linear terms come into play (for 
that we can assume a non-zero fcy). The linear terms in 
(|4l5p tend to smear or break the initial perturbation into 
wave packets propagating in opposite directions along the 
local magnetic field, such that n ~ |V^| inside those 
packets. Thus density tends to get in equipartition with 
the magnetic field and to concentrate in 2D structures as 
well. We thus expect that as a result of nonlinear stri- 
ation and linear propagation, both the density and the 
magnetic fluctuations get organized in highly intermit- 
tent, two-dimensional structures or sheets, elongated in 
z direction. This is indeed consistent with our numerical 
observations presented in Fig (2}. 

We therefore assume that essential nonlinear interac- 
tion and energy cascade take place at such 2D structures. 
Following a standard procedure [e.g.. [29|. consider tur- 
bulent fluctuations of field-perpendicular size A. Since 
such fluctuations cover 2D sheets, they occupy the vol- 
ume fraction p\ oc A. The energy density of such fluctua- 
tions therefore scales as E\ cx n\p\. The energy cascade 
time is estimated as before [cf. discussion preceding J9j] , 
t ~ 1/oj ~ l\ ~ X 3 /ipx ~ A 2 /nA, and the condition of 
constant energy flux reads E\/t — const, which gives 
n\ cx A 1 / 3 . The scaling of the energy is then E\ oc A 5 / 3 , 
and the Fourier energy spectrum scales as 

E(k ± )dk ± (xkl s/3 dk ± . (12) 

This spectrum of kinetic Alfven turbulence is in excellent 
agreement with the numerical observation in Fig. |T]), 
and it is the main result of our work. It provides a 
plausible explanation for the solar wind measurements, 
e -g- 0i EH]. Balancing the linear wave frequency 
uj ~ 1/(^A) with the inverse nonlinear interaction time 
1/t ~ n\/\ 2 ~ A -5 ^ 3 we further derive the anisotropy 
of the turbulent fluctuations with respect to the (local) 
large-scale magnetic field: I ~ A 2 / 3 . If we formally in- 
troduce the local field-parallel wave number as fc|| ~ l/l, 
then the "field-parallel energy spectrum" corresponding 

to (fl"2"|) is E{k\\)dk\\ oc fen 7 ^ 2 dk\\. This is also consistent 
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FIG. 2: Density fluctuations (top) and amplitude of magnetic- 
field fluctuations (bottom) in a field-perpendicular cross sec- 
tion of a simulation domain. The large-scale harmonics with 
k < 27r/12 have been filtered out; the plot thus represents 
fluctuations in the inertial interval, see Fig. (0. The plot 
suggests that both density and magnetic fluctuations are con- 
centrated along two-dimensional structures. Corresponding 
field-parallel cross sections (not shown here) are consistent 
with this picture. 



with the solar wind measurements [1J]. As for the 
spectrum of the electric field, it is a factor of k\ flatter 

than (HU), E E oc kj_ 2/3 H]. 

Discussion — We have proposed a model for kinetic 
Alfven turbulence below the dispersion scale (ion- 
acoustic scale). Based on numerical simulations of the 
fluid equations ([lOllip and on analytic modeling we pro- 
pose that the energy spectrum of such turbulence scales 
as , meaning that both magnetic and density fluc- 

tuations should have the same Fourier spectrum. This 



result is also consistent with in situ observations of the 
sub-proton solar wind fluctuations, and with the results 
of the gyro-kinetic simulations where the spectra close to 
kl 2 - s are observed M, [T^j • 

Our model is complementary to the previously pro- 
posed explanations invoking Landau damping, the pres- 
ence of weak kinetic Alfven turbulence, effects of wave- 
particle scattering, etc. [e.g..ll6l ll7| . These explanations 
are interesting and the effects they point out may in- 
deed affect a turbulent cascade. The difference of our ap- 
proach (as compared to the gyrokinetic numerical stud- 
ies, for example) is that it does not take into account 
Landau damping and it allows us to drive turbulence in 
a strongly coupled state. With those "spoilers" removed, 
the observed turbulent spectrum fc^ 8 ^ 3 is expected to 
arise from the nonlinear interaction, similar to the Kol- 
mogorov spectrum of hydrodynamic turbulence. This 
also implies that dissipative effects may be less important 
at the subproton scales than was previously thought. 

Our explanation points out to an interesting property 
of the kinetic- Alfven nonlinear dynamics, which tend to 
concentrate magnetic and density fluctuations at two- 
dimensional structures. This leads to strong spatio- 
temporal intermittency in the field distributions. In this 
work we have studied only the second-order statistics of 
the fluctuating fields (expressed through energy spectra) ; 
we plan to present more detailed discussion of the sub- 
proton turbulence elsewhere. 

Finally, we note that a system formally similar to 
our system of equations (ITU1) . (fTTj) also appears in the 
limit of strong guide magnetic field in the so-called elec- 
tron MHD, where the ions are assumed to be immobile. 
The wave modes in this case correspond to the so-called 
whistler waves and belong to a different branch of plasma 
dispersion relations (this branch continues into a magne- 
tohydrodynamic mode or compressional Alfven mode in 
the limit of low k, while the kinetic Alfven mode contin- 
ues into the shear Alfven mode, [e.g., H3|). One there- 
fore expects that our simulations of turbulence based 
on the system (jlOllip should be relevant for electron 
MHD in the limit of strong guide field. Previous simu- 
lations of electron MHD were conducted mostly for two- 
dimensional cases and/or for relatively weak guide fields 
and/or decaying cases, and the spectrum close to kj 7 ^ 3 
was observed [e.g., 0, d, H3, [H, [3l|. Our results suggest 
that in the limit of strong guide field, the spectrum of 
driven strong electron-MHD turbulence may be modified 
as well. 
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